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On the basis of density functional - pseudopotential calculations, we study structural, electronic, 
magnetic, and mechanical properties of the hypothetical CaC ionic compound in the rock-salt (RS), 
B2, zinc-blende (ZB), wurtzite (WZ), NiAs (NA), and anti-NiAs (NA*) structures. The results show 
that RS-CaC is the most stable system at equilibrium while applying negative hydrostatic pressures 
may stabilize the half-metallic WZ and ZB structures. The ferromagnetic equilibrium state observed 
in the RS, ZB, NA and WZ structures of CaC is attributed to the sharp partially filled p band of the 
carbon atom. It is argued that the ionic interaction increases the sharpness of the p band and hence 
enhances ferromagnetism while the covalent interaction increases the band dispersion and weakens 
magnetism. We investigate various properties of the exchange interaction in the ferromagnetic 
CaC structures. It is observed that the interatomic exchange interaction in these systems have 
consistent behavior with the spin splittings of the bond points. Comparing the structural properties 
in the ferromagnetic and nonmagnetic states, indicates a weak magneto-structural coupling in CaC. 
Applying non-local corrections to the exchange functional enhance the exchange splitting and hence 
give rise to a half-metallic electronic structure for the WZ, RS, and NA structures of CaC. In 
the framework of density functional perturbation theory, the phonon spectra of these systems are 
investigated and the observed dynamical instabilities in the NA* and B2 structures are attributed 
to the tendency of the carbon atoms toward dimerization. Various mechanical properties of the 
dynamically stable structures of CaC are determined from their phonon spectra. 



I. INTRODUCTION 

Recent observation of p magnetism in some materials 
has attracted considerable research interest as a novel 
magnetic phenomenon different from conventional mag- 
netic phenomena originating from a partially filled d or f 
shell. ^ Using first-principle calculations, in 2004, Kusak- 
abe et al^ predicted a new class of the p magnetic sys- 
tems in some hypothetical binary ionic compounds with 
no transition metals including CaP, CaAs and CaSb. 
The additional advantage of these systems is that in 
the metastable zinc-blende (ZB) phase, these compounds 
show half-metallic behavior which is attractive in the 
Spintronics community.- In a half-metal, only one spin 
channel contributes to the states around the Fermi level. 

Later on, Sieberer et al.^ predicted half-metallicity in 
the hypothetical ZB structure of several other 11"^- V"*, 
II'^-IV'^, and I'^-V'^ compounds such as CaN, CaC, and 
LiP. In all of these ionic binary compounds, the par- 
tially occupied sharp p band of the anion atom enhances 
the Stoner exchange interaction and thus give rises to 
a ferromagnetic (FM) ground state. This mechanism is 
schematically sketched in Fig. [TJ As it is qualitatively 
argued in the figure, if exchange splitting is larger than 
band width, a half-metallic gap appears in the majority 
channel. The flatness of the p band is controlled by the 
cell volume, system ionicity, and slight p-d hybridization 
around the Fermi leveL^ 

Two major open issues in the novel p magnetic com- 
pounds are their stable structure and properties of the 
exchange interaction in these systems. First-principles 
investigations may provide very helpful and accurate in- 
formation for understanding these open issues. A point 
worth to mention is that in conventional ferromagnets. 



as the highly correlated d or f electrons have the key 
role, more accurate many body calculations are usually 
needed to describe the system while in the p magnets less 
expensive exchange-correlation calculations are expected 
to be sufficient. 

In this paper, we focus on the hypothetical CaC com- 
pound as a potential p magnetic system and apply den- 
sity functional calculations to investigate its structural 
and magnetic properties in six different structures: cubic 
rock-salt (RS), cesium chloride (B2) and ZB and hexag- 
onal wurtzite (WZ), NiAs (NA) and anti-NiAs (NA*) 
structures. The NA* structure is made by permuting Ni 
and As in the NA structure. The natural compound of 
Ca and C is CaC2 that is a nonmagnetic insulator and 
crystallizes in the rock-salt structure in which C2~ dimers 
and Ca^"*" ions occupy the 4a and 4b wyckoff positions;^ 
Understanding the p ferromagnetism in the hypotheti- 
cal CaC compound may provide new opportunities for 
developing ferromagnetic carbon nanostructures. 
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FIG. 1: (Color online) Schematic spin resolved density of 
states (DOS) of a p magnetic system in the non-magnetic 
(NM), ferromagnetic (FM) and half-metallic (HM) states. 
The upward and downward plots show the majority and mi- 
nority states and the dashed lines show the Fermi energy. 
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TABLE I: Reciprocal space grids applied to the PBE, phonon 
and HSE calculations. 





PBE 


Phonon 


HSE 


RS 


10 X 10 X 10 


4x4x4 


5x5x5 


ZB 


10 X 10 X 10 


4x4x4 


5x5x5 


wz 


12 X 12 X 10 


5x5x4 


5x5x5 


NA 


15 X 15 X 9 


5x5x3 


5x5x3 


B2 


20 X 20 X 20 


6x6x6 




NA* 


16 X 16 X 16 


5x5x5 





After reviewing our computational methods in the next 
section, the structural and magnetic properties of all se- 
lected CaC structures will be discussed in section III. 
Then various properties of the exchange interaction in 
CaC compound will be analyzed in section IV. Next, the 
phonon spectra of the systems will be presented to in- 
vestigate the dynamical stability and elastic properties 
of the CaC structures. Our conclusions are presented in 
the last section. 



II. COMPUTATIONAL METHODS 

The geometry optimizations and electronic structure 
calculations are carried out by using the PWscf code 
of the QuANTUM-ESPRESSO package The gener- 
alized gradient approximation (GGA) in the scheme 
of Perdew, Burke and Ernzerhof (PBE)i and norm- 
conserving pseudopotentialSfS' by considering 10 valence 
electrons for Ca and 4 valence electrons for C, were used 
throughout our calculations. The Kohn-Sham single par- 
ticle wave functions were expanded in plane wave basis 
set up to energy cutoff of 80 Ry and the Fourier expan- 
sion of electron density were cut at 320 Ry. 

The Brillouin zone integrations were performed on the 
F-centered symmetry-reduced Monkhorst-Pack^ k-point 
meshes (table H]), using the Methfessel-PaxtoniS smear- 
ing method with a broadening parameter of 1 mRy. All 
computational parameters were optimized to achieve to- 
tal energy accuracy of about 0.1 mRy /formula unit (fu). 
The electronic structure calculations were performed in 
the scalar relativistic limit, ignoring the relativistic spin- 
orbit interaction which is expected to be small in the light 
carbon and calcium atoms. The phonon frequencies are 
calculated using the density functional linear response 
method for metallic systemsiii Table|T]lists the reciprocal 
space grids used for the expensive phonon calculations. 

We employ topological analysis of the total as well as 
spin resolved electron charge density,^- for accurate de- 
scription of bonding and magnetism in the systems. In 
this scheme, bond points are defined as the saddle points 
of electron density between neighboring nuclei. At the 
bond saddle points, electron density displays a minimum 
in the bond direction and two maxima in the perpendic- 
ular directions. 
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FIG. 2: (Color online) Calculated cohesive energy and mag- 
netic moment (per formula unit) of CaC in the six studied 
structures as a function of the primitive volume. The total 
energy diagrams are plotted in limited ranges around the equi- 
librium volumes while a broader range is used for plotting the 
diagrams of the magnetic moments. 



III. STRUCTURAL AND MAGNETIC 
PROPERTIES 

In order to find the structural properties of CaC in 
the selected structures, we computed the primitive cell 
energy of all structures as a function of the primitive 
volume, in the ferromagnetic and nonmagnetic phases. 
The B2 and NA* structures of CaC were found to be 
nonmagnetic around the equilibrium volume, while other 
systems (ZB, WZ, NA, and RS) stabilized in a ferromag- 
netic ground state. The total energy and magnetic mo- 
ment of all systems in their stable magnetic states are 
plotted as a function of volume in Fig. [2] It is observed 
that the nonmagnetic NA* and B2 structures have the 
lowest energies, while the ferromagnetic ZB phase is the 
highest energy metastable structure among the studied 
systems. The equilibrium magnetic moments of WZ, RS, 
and NA ferromagnetic structures are slightly less than 
2 [IbIJu while the ferromagnetic ZB-CaC has an inte- 
ger moment of 2 hb/ fu. The integer value of the total 
magnetic moment in the ZB structure evidences the half- 
metallic behavior of this high energy metastable struc- 
ture, while other ferromagnetic structures require some 
negative pressures (volume expansion) to achieve a half- 
metallic electronic structure. 

In order to address the relative stability of the studied 
systems at elevated pressures, we calculated and plotted 
in Fig. [31 the enthalpy free energy [H] of the systems: 
H — E — pV , where E is the converged primitive cell en- 
ergy, V is the primitive volume, and p is the hydrostatic 
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FIG. 3: (Color online) Calculated enthalpy free energy of CaC 
in the six studied structures as a function of the hydrostatic 
pressure and the resulting pressure phase diagram. 



pressure calculated from the slope of the energy- volume 
curves. The resulting simple pressure phase diagram, 
plotted in the same figure, indicates thermodynamic sta- 
bility of the NA* structure at the pressures above about 
—5.5 GPa and RS is the second stable structure in this 
region. Because of the dynamical instability of the NA* 
and B2 structures, will be shown in the section "Dy- 
namical properties" , in practice, the RS structure is ex- 
pected to be the most stable phase at the pressures above 
—5.5 GPa. Below this transition pressure, the WZ struc- 
ture becomes more stable. The pressure induced RS-WZ 
structural transition seems to be a first order transition, 
because the enthalpy of these two structures does not 
connect smoothly at the transition pressure. Inspecting 
the calculated magnetic moments (Fig. [5]) indicates that 
the RS-WZ structural transition is accompanied by a nor- 
mal ferromagnet to half-metal phase transition. The ZB 
structure finds thermodynamic stability after a further 
transition at a pressure of about —6.5 GPa. Although, 
the obtained phase diagram shows the stability of the 
half-metallic systems in the negative pressures, advanced 
epitaxial growth techniques make it feasible to synthe- 
size the stretched thin films of these systems on proper 
substrates. 

The calculated equilibrium properties determined by 
fitting the Murnaghan equation of statei2, to the calcu- 
lated total energy- volume data are presented in Fig. SI 
It is seen that in the sequence of ZB, WZ, NA, RS, B2, 
and NA* structures, the cohesive energy, equilibrium vol- 
ume, and compressibility (inverse of the bulk modulus) 
are decreasing. These behaviors indicate enhancement 
of the bond strength and bond stiffness in this sequence. 
For better understanding, the ionization parameter / and 
the total electron density at bond points pB (per formula 
unit) are also computed and plotted in Fig. S) The ion- 
ization parameter, defined as the transferred charge from 
cation to anion normalized by the cation valence charge, 
is a measure of the ionic interaction in the systems. On 
the other hand, is the sum of electron densities at 
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FIG. 4: (Color online) Calculated equilibrium properties of 
CaC in the selected structures. V: equilibrium primitive vol- 
umes, B: bulk modulus, Ec'- cohesive energy, Em- magneti- 
zation energy defined as the difference between the ferromag- 
netic and nonmagnetic cohesive energies, /: ionization, pa: 
total electron density at bond points. 



all topological bond points (per formula unit) and thus 
speculates the strength of the covalent bonding in the 
systems. We observe that the enhancement of the bond 
strength and stiffness in the studies structures is accom- 
panied by decreasing the ionic and increasing the cova- 
lent interactions in the systems. It is known that the di- 
rectional covalent bonds form solids with high strength, 
high stiffness, and high melting point while ionic inter- 
action typically give rises to intermediate strength and 
extremely brittle solids. In the next sections, we qual- 
itatively compare the brittleness of the systems. The 
calculated results (Fig. 2]) indicate that the nonmagnetic 
B2 and NA* structures exhibit the highest covalency and 
the lowest ionization among the studied systems. More- 
over, we observe that the more ionized structures (except 
NA) have larger absolute magnetization energies. These 
observations are evidences for the ionic origin of ferro- 
magnetism in CaC. 

The equilibrium lattice parameters, Fermi level spin 
polarizations, bond lengths, and atomic magnetic mo- 
ments in the six studied lattices of CaC are listed in ta- 
ble ini After topological analysis of charge densities, we 
identified that, in addition to the ionic carbon-calcium 
bonding in all systems, there are some covalent carbon- 
carbon bonds in the nonmagnetic NA* and B2 structures. 
Therefore, equilibrium properties of the C-C bond, in 
addition to C-Ca, are also listed in table HIl In order to 
understand these results, we note that while the C-Ca in- 
teraction is mainly ionic, the C-C bonding is totally cova- 
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TABLE II: Computed equilibrium lattice parameters a and c, atomic and total magnetic moments, and C-C and C-Ca bond 
properties in the selected structures of CaC; P{%): total spin polarization at the Fermi level, d(A): equilibrium distance, Ni,: 
number of the topological bonds per formula unit, and pb (e/A^): bond electron density. The WZ structure involves slightly 
different vertical and in-plane bonds and hence, in this case, the average bond properties are reported. 
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3.64 







2.65 


6 
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WZ 


4.38 


5.30 


90 
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0.21 


1.99 
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0.20 


ZB 


5.75 




100 


1.86 


0.19 


2.00 


4.07 







2.49 


4 


0.24 



lent. Observation of the topological covalent C-C bonds 
in the nonmagnetie B2 and NA* structures explains the 
highest covalency of these systems (Fig. U]). Moreover, 
the slight increasing of the magnetic moment in RS, NA, 
WZ, and ZB series is generally accompanied by increas- 
ing (decreasing) the C-C (C-Ca) distance. Decreasing 
the C-Ca distance intensifies the ionic attraction in the 
systems and hence enhances the magnetic moment via 
increasing the charge transferred from Ca to the flat C p 
orbital (Fig.|31). In contrary, the reduction of the C-C dis- 
tance enhances the covalent interaction in the system and 
consequently increases the valence band width and then 
decreases the magnetic moment. These findings, further 
evidence that ferromagnetic behavior of CaC roots into 
the ionic C-Ca interaction. 

All magnetic structures show a magnetic moment of 
about 2 /isZ/u with no topological covalent C-C bond. 
In the ZB structure, the magnetic moment is exactly 
equal to 2 which is due to the half-metallic be- 

havior of this system argued by the calculated perfect 
Fermi level spin polarization (table HI)). We observe that 
the magnetic moments are mainly carried by the C atoms 
and the Ca atoms have a small moment parallel to the 
C moment. The dominant magnetic contribution of car- 
bon comes from the strong ionic C-Ca interaction which 
highly evacuates the valence s shell of calcium. The 
empty spin down C p orbital give rises to more charge 
transferring in the minority channel and leaving a posi- 
tive moment on the Ca atom. All FM structures of CaC 
have a perfect or nearly semiconducting majority chan- 
nel and the Fermi level spin polarizations comes form the 
minority electrons. Hence it seems that these systems are 
capable of conducting spin currents antiparallel to their 
magnetic moment. 

It was mentioned that all systems have C-Ca bonds 
while the topological C-C bonds appear only in the non- 
magnetic NA* and B2 structures (table HIl. In the NA* 
structure, strong covalent bonding between C atoms sig- 
nificantly decrease the c/a ratio while in the NA lattice 
the mainly ionic C-Ca interaction increases this ratio. A 
visible trend is that systems with higher number of bonds 
per atom (Nb) have generally smaller electron densities 



at bond points (pb) and longer bond lengths. It may be 
explained in this way that in these systems the valence 
electrons are divided into more bonds, hence each bond 
receives less electrons and consequently gets weaker and 
longer. However, overall bonding in the highly coordi- 
nated systems is stronger, because these systems have 
generally lower cohesive energy and higher total electron 
density at bond points (Fig. |4]). 

The heat of formation (AHf) of bulk CaC which mea- 
sures the required energy to decompose this compound 
into graphite and bulk fee Ca, is calculated as follows: 

ATT TpCaC I TpC 1 TpCa \ 

ZlHf - h^^f^ - [l^coh + ^coh) 

where the right hand side parameters are the cohesive 
energies of bulk CaC, graphite, and bulk fee calcium. 
The calculated formation energies of CaC in the studied 
structures are 1.83 eV (NA*), 2.1 eV (B2), 2.1 eV (RS), 
2.3 eV (NA), 2.4 eV (WZ) and 2.6 eV (ZB). These pos- 
itive formation energies indicate that synthesis of these 
systems requires endothermic processes. 



IV. EXCHANGE INTERACTION IN CaC 

As it was mentioned, an important issue in the novel 
field of p magnetism, is understanding the properties of 
the exchange interaction in these systems. In this regard, 
we investigate various spin dependent parameters in CaC 
ferromagnetic structures (ZB, WZ, RS, NA) to discuss 
some features of the exchange interaction in this system. 



A. Exchange stiffness 

The first calculated parameter is the energy difference 
between the ferromagnetic and antiferromagnetic (AF) 
spin states of CaC in the selected structures. Applying a 
simple effective Heisenberg model to the exchange inter- 
action between the nearest neighbor carbon atoms with 
magnetic moments Mi and M2: E = J/2(Mi • M2) en- 
ables us to estimate the interatomic exchange stiffness 
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FIG. 5: Calculated exchange stiffness J (meV), bond spin 
splitting Abond (mA), and exchange splitting of the carbon 
p states ADOS (eV) in the ferromagnetic CaC structures. 



from the energy difFerence of the FM and AF states: 
J = {Eaf — Epm)/\M\'^. The estimated J parameters 
may be used for qualitative comparison of the Curie tem- 
peratures of different systems. In the hexagonal NA and 
WZ structures, the AF spin ordering were applied in the 
[0001] direction while in the cubic RS and ZB structures 
we calculated and compared the minimized energy of the 
type I (AF ordering in the [001] direction) and type II 
(AF ordering in the [111] direction) antiferromagnetism 
and found that both structures prefer AF-type II. Em- 
ploying the obtained antiferromagnetic energies, the in- 
teratomic exchange stiffness between the nearest C atoms 
in the different CaC structures was estimated to be about 
24 meV (RS), 42 meV (NA), 52 meV (WZ), and 49 meV 
(ZB). Hence, it is speculated that WZ and ZB have the 
highest Curie temperatures among these four structures. 
The obtained J values, plotted in Fig. [SI are substan- 
tially larger than the exchange stiffness of some strong 
ferromagnets such as Co2MnSi Heusler alloyJ^ Then, 
RS-, WZ-, ZB-, and NA-CaC are expected to be high 
Curie temperature ferromagnets and thus promising for 
practical applications. The theoretical Curie tempera- 
ture of ZB-CaC, in the mean field and random phase 
approximations, is reported to be about 1051 and 735 K, 
respectively!^ 



B. Exchange splittings 

The spin resolved partial electronic density of states 
(DOS) of the valence states of CaC in its ferromagnetic 
structures are displayed in Fig. [6l It is seen that the car- 
bon p states have the dominant contribution to the va- 
lence band with slight hybridization with the calcium d 
states. Therefore, magnetization mainly comes from the 
carbon p electrons (Table. [Tl]). The valance bands in the 
ZB structure are considerably sharper, compared with 
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FIG. 6: (Color online) Calculated partial DOS of the carbon p 
(solid line) and calcium d (shaded areas) orbitals in the ferro- 
magnetic CaC structures. The positive and negative densities 
correspond to the majority and minority states, respectively. 
The Fermi energies are set to zero. 



the other three structures. This observation explains the 
half-metallic behavior of ZB-CaC in equilibrium volume. 
By using the sharp peaks of the partial DOS in the ma- 
jority and minority channels, we estimated the exchange 
splitting of the carbon p states in all ferromagnetic struc- 
tures and plotted the obtained data in Fig. [SJ It is clearly 
visible that the exchange splitting of the partial p DOS 
is not correlated with the exchange stiffness parameter. 
The reason is that the interatomic exchange coupling 
constant J only measures the exchange interaction be- 
tween the neighboring carbon valence electrons while the 
exchange splitting of the carbon p states originates from 
the full exchange interaction between all valence states. 
In other words, the on-site exchange interaction plays an 
important role in the exchange splitting of the carbon 
p states while it has small effect on the interatomic ex- 
change coupling. An evidence for this fact is that, when 
the interatomic distances increase to large values, the on- 
site exchange splittings are usually enhanced while the 
interatomic exchange couplings vanish. 

Since the bond points are rather far from the nuclei, 
their spin dependent properties are expected to be less 
sensitive to the on-site exchange interaction and thus bet- 
ter correlated with the interatomic exchange coupling. 
In order to verify this statement, we applied topological 
analysis to the spin resolved electron densities to find the 
majority and minority bond positions in the ferromag- 
netic structures. The distance between the spin resolved 
topological bond points is defined as the bond spin split- 
ting and plotted in Fig. [S] As it was stated, we observe a 
clear correlation between this distance and the exchange 
stiffness. Therefore, the bond spin splitting may be as- 
sumed as the real space representation of the interatomic 
exchange interaction in ferromagnetic materials. 
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TABLE III: Calculated cohesive energy (Ec), equilibrium volume (Vo), bulk modulus (B), ionization (/), and bond electron 
density (pb) of the ferromagnetic structures of CaC and the ZB structure of MnAs and CrAs in the nonmagnetic (NM) and 
ferromagnetic (FM) phases. 
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B (GPa) 




1 (/o) 


Pb (e/A ) 




Lattice 


FM 


NM 


FM 


NM 


FM 


NM 


FM 


NM 


FM NM 


CaC 


RS 


-7.53 


-7.21 


37.1 


36.0 


60 


60 


67 


65 


0.18 0.18 




NA 


-7.31 


-7.16 


37.2 


35.2 


57 


65 


67 


66 


0.18 0.19 




WZ 


-7.24 


-6.85 


44.1 


44.3 


48 


47 


68 


68 


0.20 0.20 




ZB 


-7.00 


-6.38 


47.5 


46.5 


42 


41 


68 


68 


0.24 0.25 


MnAs 


ZB 


-6.28 


-5.39 


46.1 


37.3 


50 


105 






0.43 0.59 


CrAs 


ZB 


-6.35 


-5.30 


45.3 


39.7 


67 


96 






0.46 0.55 



C. Magneto-structural coupling 

In order to study the interplay between magnetism 
and bonding, we calculated the structural properties of 
the ferromagnetic structures of CaC in the nonmagnetic 
state. The results, presented in table Hill point to an in- 
teresting property of this class of magnetic materials; the 
equilibrium volume and bulk modulus are nearly inde- 
pendent of the magnetic state. It is then concluded that 
the exchange interaction in the CaC structures weakly af- 
fects the structural parameters of the system. For more 
understanding, we calculated the equilibrium volume and 
bulk modulus of ZB-MnAs and ZB-CrAs in the ferromag- 
netic and nonmagnetic states (table IIIip and observed 
that the nonmagnetic systems have considerably higher 
bulk modulus and lower equilibrium volume. Since mag- 
netism in MnAs and CrAs originates from the partially 
filled 3d states, it is concluded that the interplay between 
bonding and magnetism in the p magnetic systems, in 
contrast to the d magnetic systems, is weak. 

For understanding the magneto-structural coupling, 
one should take into account the way the exchange in- 
teraction infiucnces the spatial distribution of electrons 
in a magnetic system. The exchange interaction, in gen- 
eral, transfers some electrons from the occupied minor- 
ity states into the unoccupied majority states. Because 
of the symmetry splitting of the d orbital in all crys- 
tal structures, the occupied and unoccupied d valence 
states have always different orbital character and spa- 
tial distribution. For example, in the cubic crystal field, 
the d orbital splits into the triply degenerate t2 and 
the doubly degenerate e states with different spatial dis- 
tribution. Among these two sub-orbitals, the one that 
contributes more in the interatomic bonds has lower en- 
ergy and hence is more occupied while the higher energy 
and less occupied sub-orbital distributes in the out-of- 
bond regions of the lattice. Therefore, the exchange in- 
duced charge transfer in the d magnetic systems usually 
weakens the bonds by driving the minority electrons of 
the interatomic bonds toward the out-of-bond regions of 
the lattice where the majority unoccupied states are dis- 
tributed. As a result of that, the exchange interaction 



in these systems acts as an effective repulsive force and 
hence increases the interatomic distances and compress- 
ibility of the systems. 

On the other hand, the p orbital is less influences by 
the crystal symmetries. We focus on the cubic symmetry, 
in which the spin resolved p orbitals are triply degenerate 
and hence the occupied and unoccupied p valence states 
have similar spatial distributions. As a results of that, 
in the cubic p magnetic systems, the exchange induced 
charge transfer preserves the spatial distribution of elec- 
trons and consequently the nonmagnetic and magnetic 
phases of these systems have close interatomic distances 
and compressibilities. In other words, the exchange inter- 
action in the p magnetic materials, mainly influences the 
local spin distribution of electrons without signiflcant ef- 
fect on the total electron density. In order to verify this 
statement, we calculated the ionization and bond elec- 
tron density of RS-, NA-, WZ-, and ZB-CaC as well as 
ZB-MnAs and ZB-CrAs in the nonmagnetic and ferro- 
magnetic states (table Uni. It is clearly seen that in the 
CaC structures, nonmagnetic and ferromagnetic systems, 
have very similar ionizations and bond electron densities, 
while in the d ferromagnet CrAs and MnAs, the exchange 
interaction effectively reduces the bond electron densities 
of the ferromagnetic states, compared with the nonmag- 
netic states. 

The splitting of the p orbital into Pz and Pxy, in the 
hexagonal symmetry, may enhance magneto-structural 
coupling in the hexagonal p magnetic systems. In the 
hexagonal NA-CaC, Pz has lower energy and more con- 
tribution to bonding. Hence, the exchange interaction 
transfers some minority Pz electrons into the empty ma- 
jority Pxy states and consequently enhances bonding in 
the xy plane. As a result of that, NA-CaC in the fer- 
romagnetic state, compared with the nonmagnetic state, 
has smaller in-plane and larger vertical lattice parame- 
ter, giving rise to a slightly larger equilibrium volume and 
compressibility for the ferromagnetic state. On the other 
hand, in the hexagonal WZ-CaC, the exchange induced 
in-plane bond enhancement and vertical bond weaken- 
ing compensate each other and hence the FM and NM 
phases of this system have close equilibrium volume and 
bulk modulus. 
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TABLE IV: Calculated equilibrium volume (Vb), bulk modulus (B), cohesive energy (Ec), magnetic moment (/i), spin flip gap 
(sgap), and exchange splitting of DOS (ADOS) in ferromagnetic structures of CaC within the HSE03 and PBE functionals. 





Vo (A^) 


B (GP 


a) 


Ec (eV) 




) 


sgap (eV) 


ADOS (eV) 




HSE 


PBE 


HSE 


PBE 


HSE 


PBE 


HSE 


PBE 


HSE 


PBE 


HSE PBE 


RS 


37.2 


37.1 


62 


60 


-7.06 


-7.53 


2.00 


1.97 


0.8 





3.35 1.65 


NA 


37.2 


37.2 


59 


57 


-6.84 


-7.31 


2.00 


1.98 


0.8 





2.90 1.61 


WZ 


44.1 


44.1 


49 


48 


-6.78 


-7.24 


2.00 


1.99 


0.8 





3.00 1.70 


ZB 


47.6 


47.5 


42 


42 


-6.51 


-7.00 


2.00 


2.00 


1.8 


0.9 


2.95 1.65 



D. Non-local exchange effects 

The localized nature of the conventional exchange- 
correlation functionals give rises to the typical under- 
estimation of the semiconductor band gaps, compared 
with the experimental dataJ^ Therefore, because of the 
semiconducting behavior of half-metals in one spin chan- 
nel, the non-local exchange effects may significantly infiu- 
ence the half-metallic properties. The calculated Fermi 
level spin polarizations (table |ll]) and density of states 
(Fig. ini) within the localized PBE functional, indicate 
half-metallic behavior of ZB-CaC at its equilibrium vol- 
ume and nearly half-metallic character of the other ferro- 
magnetic structures of CaC. Thus applying a more accu- 
rate non-local exchange-correlation functional may sig- 
nificantly improve the observed half-metallic and nearly 
half-metallic properties. 

We use the screened HSE03 hybrid functional to apply 
the non-local exchange corrections to our calculations.— 
In this scheme the short rang (SR) part of the Hartree- 
Fock exchange is mixed with the short range PBE ex- 
change while the long range (LR) part of the exchange 
and whole correlation is taken from the PBE functional: 

jT^HSE 1 j^HF.SR , 3 j-,PBE,SR , J7^PBE,LR , t^PBE 
^xc - -^^x + -^^x + ^x + 

The exclusion of the long range Hartree-Fock exchange 
from the screened HSE03 hybrid functional efficiently im- 
proves its performance in metallic systems. Moreover, 
this feature substantially reduces the required reciprocal 
space grids for performing the expensive non-local HSE 
calculations (table H]). 

The calculated structural and magnetic properties of 
CaC in its four ferromagnetic structures within the HSE 
hybrid functional are presented in table IIVI along with 
the corresponding PBE results for comparison. The pa- 
rameter spin flip gap in a half-metal denotes the energy 
distance between the valence band maximum in the semi- 
conducting spin channel and the Fermi level and hence 
measures the minimum energy required to flip spin of an 
electron in the semiconducting spin channel. This pa- 
rameter may be used to qualify half-metals, because if it 
vanishes, the half-metallicity breaks. As it was expected, 
we observe that the non-local exchange effects efficiently 
improve spin flip gap and exchange splitting of DOS and 
hence give rises to the half-metallic behavior of all CaC 
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FIG. 7: Calculated phonon band structure of CaC in the B2, 
RS, ZB, NA*, NA, and WZ structures. The bands are plotted 
in comparable scales. 



ferromagnetic structures within HSE. 

Comparing the calculated equilibrium volumes and 
bulk moduli within the HSE and PBE functionals pro- 
vides further evidence for the weakness of the magneto- 
structural coupling in the p magnetic systems. It is 
clearly seen that the non-local exchange interactions, in 
spite of the signiflcant influence on the magnetic param- 
eters, have negligible effect on the equilibrium volumes 
and bulk moduli. The HSE cohesive energies are found 
to be about 0.4 eV higher than the PBE values (table lTvT) . 
It may be attributed to the absence of the self interaction 
in the Hartree-Fock exchange energy. 

V. DYNAMICAL PROPERTIES 

The dynamical properties of CaC in the cubic RS, ZB, 
and B2 and hexagonal WZ, NA, and NA* structures are 
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FIG. 8: (Color online) Total phonon DOS (solid line) and 
carbon contribution to the phonon DOS (shaded area) in the 
B2 and NA* structures of CaC. The phonon DOS of the B2 
structure in the negative (imaginary) frequencies is enlarged 
50 times (inset) to be more visible. 



studied by calculating the phonon spectra of these sys- 
tems in the framework of density functional perturba- 
tions theoryjii The obtained phonon band structures in 
the equilibrium ground states are displayed in Fig. [T] 
The negative values in these plots correspond to the 
imaginary phonon frequencies. It is observed that the 
B2 and NA* structures which are the lowest energy sys- 
tems among the studied structures, have some imaginary 
phonon frequencies and hence are dynamically unstable. 
The dynamical instability of the B2 structure occurs in 
some transverse acoustic phonon modes in a small por- 
tion of the BZ in the [110] direction while NA*-CaC has 
a broad instability coming from an optical phonon band 
in the whole BZ and a longitudinal acoustic band in the 
[0001] direction. 

In order to understand the origin of the dynamical in- 
stabilities, we calculated the atom resolved phonon den- 
sity of states of the B2 and NA* structures (Fig. [8|). 
It is clearly seen that the imaginary phonon frequen- 
cies mainly come form carbon oscillations. This obser- 
vation indicates that the dynamical instability of the B2 
and NA* structures is due to the tendency of the car- 
bon atoms for reconfiguration. Existence of the carbon 
dimers in the natural CaC2 compound,^ provides an ex- 
perimental evidence for the tendency of carbons in CaC 
to get close together and dimerize. This statement, is ver- 
ified by the topological investigation of the electron den- 
sity of the B2 and NA* structures. The results, partially 
displayed in Fig. |9l indicates formation of the topologi- 
cal bonds between the carbon atoms in these two struc- 
tures. As it is visible in the figure, in the B2 structure, 
the coplanar covalent carbon-carbon bonds may give rise 
to reconfiguration (dimerization) of the carbon atoms in 
the diagonal [110] (and other equivalent) directions. It 
is consistent with the observed dynamical instability of 
B2-CaC in the [110] FW direction (Fig. [71). Similar ar- 
guments apply to the phonon instabilities of NA*-CaC. 

The phonon spectra of materials is a rich source of 
information on mechanical and thermodynamical prop- 




FIG. 9: (Color online) Left: Calculated topological bonds in 
B2-CaC, The green balls with letters b and B stand for the 
ionic Ca-C and covalent C-C bond points, respectively. Right: 
Helmholtz free energy of the ferromagnetic CaC structures as 
a function of temperature. 

erties. The vibrational based thermodynamic properties 
of a system is studied by calculating the vibrational con- 
tribution to the Helmholtz free energy (F) as follows: 

F(T) ^E-ksT j duD{u^) ln(^--^) 

where E is the total energy of the system, the integra- 
tion is over the phonon frequencies, D(uj) is the phonon 
DOS, kB is the Boltzmann constant, T is the kelvin tem- 
perature, and P = l/ksT. By inserting the calculated 
phonon DOS into the above equation, the Helmholtz free 
energy of CaC in its four ferromagnetic structures were 
calculated and presented in (Fig. |9l). The results indi- 
cate absence of any temperature induced structural phase 
transition in the ferromagnetic structures of CaC. 

Finally, we have used our first-principle phonon spec- 
tra to calculate various mechanical properties of the fer- 
romagnetic structures of CaC (table |V|) . The Debye tem- 
perature (Tb), which measures hardness of a solid, was 
calculated by using the speed of sound and density of 
atoms in the system.^ The elastic constants are cal- 
culated from the slopes of the acoustic phonon bands 
at the r point fi2, The obtained elastic constants were 
used to calculate several mechanical parameters, listed 
in table [V] The cubic and hexagonal bulk B and shear 
G moduli are calculated from the elastic constantsj^ 
and then the Young modulus {Y) is calculated by y = 
9BG/(3B -I- G). In the cubic ZB and RS lattices, we ob- 
serve that the phonon based bulk and Young moduli are 
close to the directly calculated values (table |Vl. This con- 
sistency confirms reliability and accuracy of our phonon 
frequency calculations. In the hexagonal structures, two 
schemes proposed by Reuss and Voigt,^ are applied for 
calculating the bulk and shear moduli from the elastic 
constants. In the NA structure we found that the Reuss 
scheme gives better bulk and Young moduli, compared 
with the directly calculated values, while in the hexago- 
nal WZ structure, the Voigt scheme exhibits better per- 
formance. The B/G ratio, presented in table |Vl qualita- 
tively determine the ductility of a material.— Therefore, 
among the ferromagnetic structures of CaC, the most 
ductile and brittle systems are expected to be the ZB 
and RS structures, respectively The mechanical parame- 
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TABLE V: Calculated mechanical properties of CaC in the RS, NA, WZ, and ZB structures, ro(K): Debye temperature, 
Cij(GPa): elastic constants, Bph, Yph, and Gph(GPa): bulk. Young, and shear moduli calculated from the elastic constants, B 
and F(GPa): directly calculated bulk and Young moduli from the energy-volume and energy-lattice parameter curves, and o: 
poison ratio. 



Structure 


To 


Cii 


C33 


C44 


Cee 


C12 


Cl3 


Bph 


B 


Yph 


Y 


Gph 


B/G 


a 


RS 


470 


153 




58 




13 




60 


60 


139 


124 


63 


0.95 


0.08 


NA 


466 


131 


196 


55 


72 


-14 


13 


54 


57 


138 


149 


64 


0.84 


0.05 


WZ 


459 


76 


200 


10 


13 


49 


2 


49 


48 


55 


57 


27 


1.81 


0.27 


ZB 


458 


53 




20 




36 




42 


42 


42 


44 


16 


2.63 


0.40 



ter poison ratio, calculated by cr = 1/2(1 — Y/3B), may 
be qualitatively used to compare structural stiffness of 
materials; more stiff materials have lower poison ratios. 
Hence the RS and NA structures of CaC are expected 
to be more stiff systems, compared with the WZ and ZB 
structures. 



VI. CONCLUSIONS 

In this paper we employed density functional- 
pseudopotential calculations to study structural, mag- 
netic, and mechanical properties of the zinc-blende (ZB), 
wurtzite (WZ), rock-sah (RS), B2, NiAs (NA), and anti- 
NiAs (NA*) structures of CaC compound. Considering 
the transferred charge from cation to anion as well as 
the electron density at bond points, we found that the 
covalency of the system increases along the ZB-WZ-NA- 
RS-B2-NA* sequence, while the ionicity decreases. As 
a result of that, the NA* and B2 structures exhibit the 
lowest cohesive energies and nonmagnetic ground states 
while the ZB structure is the highest energy metastable 
system with the half-metallic property. Investigation of 
the enthalpy free energy of the systems along with the ob- 
served dynamical instabilities of the NA* and B2 struc- 
tures showed that the ferromagnetic RS structure is the 
most stable system at the equilibrium conditions while 



stabilization of the half-metallic ZB and WZ structures 
requires some negative hydrostatic pressures (lattice ex- 
pansions). The dynamical instabilities of the nonmag- 
netic NA* and B2 structures were attributed to the ob- 
served tendency of the carbon atoms toward reconfigura- 
tion (possibly dimerization) . The observed p ferromag- 
netism in the ZB, RS, WZ, and NA structures of CaC 
was argued to be due to the strong ionic interaction be- 
tween calcium and carbon which give rises to a sharp par- 
tially filled p band on the carbon atom. The interatomic 
exchange stiffness of the ferromagnetic CaC structures 
were speculated to be higher than 20 meV. Because of 
the less influence of the crystal symmetries on the p or- 
bital, the magneto-structural coupling in the p magnetic 
systems was argued to be weak. The application of the 
non-local exchange interaction via the screened HSE03 
hybrid functional resulted in a perfect Fermi level spin 
polarization and half-metallic behavior for all ferromag- 
netic CaC structures. 
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